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Abstract 



00^ 

pH I We present a method devised by Jacobi to derive Lagrangians of any second-order 

differential equation: it consists in finding a Jacobi Last Multiplier. We illustrate the 
easiness and the power of Jacobi's method by applying it to the same equation studied 
by Musielak et al. with their own method [Musielak ZE, Roy D and Swift LD. Method 
to derive Lagrangian and Hamiltonian for a nonlinear dynamical system with variable 
coefficients. Chaos, Solitons & Fractals, 2008;58:894-902]. While they were able to find 
one particular Lagrangian after lengthy calculations, Jacobi Last Multiplier method yields 
0^ ' two different Lagrangians (and many others), of which one is that found by Musielak et al, 

and the other(s) is(are) quite new. 



O ; 1 Introduction 

oo 

O 

It should be well-known that the knowledge of a Jacobi Last Multiplier always yields a La- 
grangian of any second-order ordinary differential equation [9], [27j. Yet many distinguished 
p\ ' scientists seem to be unaware of this classical result. In this paper we present again the method 

j^ . of the Jacobi Last Multiplier in order to compare the easiness and the power of Jacobi's method 



with that proposed by Museliak et al |12] for the same purpose. We have already presented 
Jacobi's method in [T7]. The references in [17] and the papers |15j-|26j may give an idea of the 
many fields of applications yielded by Jacobi Last Multiplier. 

In |12j the authors searched for a Lagrangian of the following second-order ordinary differ- 
ential equation 

X + b{x)x + c{x)x = (1) 
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with b{x),c{x) arbitrary functions of the dependent variable x = x{t). After some lengthy 
calculations they found one Lagrangian. In the present paper we apply Jacobi's method to 
equation ([1]) and show that many (an infinite number of) Lagrangians can be easily derived. 

This paper is organized in the following way. In section 2, we illustrate the Jacobi last 
multiplier and its properties [S]-[S], its connection to Lie symmetries [TU], |llj . and its link 
to the Lagrangian of any second-order differential equations [9], [27]. Moreover we exemplify 
Jacobi's method by determining the Lagrangians of two equations studied by Euler ^ and 
Jacobi himself f8] for the purpose of finding their multipliers. In section 3, we apply Jacobi's 
method to equation ([1]) and determine some of its many Lagrangians. In section 4, we conclude 
with some final remarks. 

Here we employ ad hoc interactive programs [H] written in REDUCE language to calculate 
the Lie symmetry algebra of the equations we study. 



2 The method by Jacobi 

The method of the Jacobi last multiplier [6]- [9]) provides a means to determine all the solutions 
of the partial differential equation 



Af = y^^ai{xi,...,Xn) 
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or its equivalent associated Lagrange's system 
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In fact, if one knows the Jacobi last multiplier and all but one of the solutions, then the last 
solution can be obtained by a quadrature. The Jacobi last multiplier M is given by 



d{f,u;i,uJ2,---,uj, 



n-l 



where 





d{xi,X2,...,Xn) 


^J : 






r df 


df 


d{f,UJl,UJ2,- 
d{xi,X2,. 


•'^7^)-det 


dxi 
doji 

dxi 


dXn 
OUJi 

dXn 






duJn-l 


duJn-l 



dxi 



dXn 



(4) 



(5) 



and wi, . . . ,uJn-i are n — 1 solutions of ([2]) or, equivalently, first integrals of ^ independent 
of each other. This means that M is a function of the variables (xi, . . . ,x„) and depends on 
the chosen n — 1 solutions, in the sense that it varies as they vary. The essential properties of 
the Jacobi last multiplier are: 

(a) If one selects a different set of n — 1 independent solutions rji,. . . ,r]n-i of equation ([2]), 
then the corresponding last multiplier A^ is linked to M by the relationship: 



d{LL>i,...,UJn-l) 

(b) Given a non-singular transformation of variables 

T : (Xi, X2, • • • , Xnj > \^li ^2) • • • ' ■^nJj 

then the last multiplier M' of A'F = is given by: 

^d{xi,X2,...,Xn) 



M' = M 



uyx^, ^2, • • • , X„J 



where M obviously comes from the n — 1 solutions of AF = which correspond to those 
chosen for A'F = through the inverse transformation r^^. 

(c) One can prove that each multiplier M is a solution of the following linear partial differ- 
ential equation: 

E^^O; (6, 

i = l * 

viceversa every solution M of this equation is a Jacobi last multiplier. 

(d) If one knows two Jacobi last multipliers Mi and M2 of equation ([2]), then their ratio is a 
solution tj of ([2]), or, equivalently, a first integral of ([3]). Naturally the ratio may be quite 
trivial, namely a constant. Viceversa the product of a multiplier Mi times any solution 
uj yields another last multiplier M2 = Miuj. 

Since the existence of a solution/first integral is consequent upon the existence of symmetry, 
an alternate formulation in terms of symmetries was provided by Lie [llj. A clear treatment 
of the formulation in terms of solutions/first integrals and symmetries is given by Bianchi |3] . 
If we know n — 1 symmetries of ([2])/([3|), say 

n 

ri = '^^ij{xi,...,Xn)dx^, i = l,n-l, (7) 



Jacobi's last multiplier is given by M = A , provided that A 7^ 0, where 
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There is an obvious corollary to the results of Jacobi mentioned above. In the case that there 
exists a constant multiplier, the determinant is a first integral. This result is potentially very 
useful in the search for first integrals of systems of ordinary differential equations. In particular, 
if each component of the vector field of the equation of motion is missing the variable associated 
with that component, i.e., dai/dxi = 0, the last multiplier is a constant, and any other Jacobi 
Last Multiplier is a first integral. 

Another property of the Jacobi Last Multiplier is its (almost forgotten) relationship with 
the Lagrangian, L = L{t,x,x), for any second-order equation 



is [9], [27] 



X = F{t, x,x) 
M 



dx"^ 



where M = M{t, x, x) satisfies the following equation 

d rlW 

_(,o,M) + _ = 0. 

Then equation Q becomes the Euler-Lagrangian equation: 

_ A f^\ —-(] 
dt \dx J dx 
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The proof is given by taking the derivative of p2|) by x and showing that this yields pi|) . If 
one knows a Jacobi last multiplier, then L can be easily obtained by a double integration, i.e.: 



Mdx] dx + fi{t,x)x + f2{t,x), 



(13) 



where /i and /2 are functions of t and x which have to satisfy a single partial differential 
equation related to ([9]) [22]. As it was shown in [22], /i,/2 are related to the gauge function 
g = g{t,x). In fact, we may assume 
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where /s has to satisfy the mentioned partial differential equation and g is obviously arbitrary. 
In [22] it was shown that if one knows several (at least two) Lie symmetries of the second- 
order differential equation ([9]), i.e. 



r^- = Vjit, x)dt + Gjit, x)d^, j = 1, r, 
then many Jacobi Last Multipliers could be derived by means of (| 



(15) 



I.e. 
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with (n, ?TT, = l,r), and therefore many Lagrangians can be obtained by means of (I13p . 

In |22j fourteen different Lagrangianq^ were derived even for an equation as controversial 
as the damped linear harmonic oscillator, i.e.: 



X + ex + /ex = 0, 



(17) 



which about 80 years ago was thought not to be derivable from a variational principle by Bauer 
[2] and to be "a physically incomplete system" , namely in need of additional equationqj by 
Bateman pj. 



2.1 Two examples by Euler in jS] 

In ^, Jacobi found his "new multiplier" for the following class of second-order ordinary dif- 
ferential equation^ studied by Euler [4J [Sect. I, Ch. VI, §§915 ff.]: 



1 5c/7 . n dip . ^ 
x+--^x2 + -^x + S 
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(18) 



with (/?, B arbitrary functions of t and x. Indeed Jacobi derived that the multiplier of equation 
p^ is given by: 



M 



oV{t,3:) 



(19) 



as it is obvious from (jlip . Then, Jacobi presented two examples of the class of equations (|18p . 
also studied by Euler, to illustrate the use of his multiplier, namely that the knowledge of 



^The fourteen Lagrangians are independent from each other and not related by any gauge function. They are 
derived from the eight-dimensional Lie symmetry algebra which is admitted by any linear second-order ordinary 
differential equations [llj . 

^Bateman called them "a complementary sets of equations". 

^This is not Jacobi's original notation. 



one first integral and one multiplier yields integration by quadrature. Here we report those 
equations for the reader's conveniencqj- We also find the corresponding Lagrangian and search 
for Lie symmetries of Euler's equation in order to describe the link between Lie symmetries, 
Jacobi Last Multiplier, first integrals, and of course Lagrangians. In fact it easy to derive a 
Lagrangian of equation (flBj) by means by means of PU]) . i.e. 

L = ^e^(t'-)x' + f,(t,x)± + f2{t,x) (20) 

with /i, /2 functions of t and x satisfying the following equation: 

2.1.1 Example I 

The first example by Euler that Jacobi presented is the following equation: 

x'^x + xx^ + /3x - 7t = (22) 

with (3 and 7 arbitrary constants. Jacobi obtained the multiplier 

Ml = x^, (23) 



by means of (J19p and showed how to integrate equation (j22p . 

We can use Jacobi's multiplier ([23]) to derive a Lagrangian of equation ([22]) by means of 
(UnD, i.e. 

Li = -xV + /i(t,x)x + /2(t,x) (24) 

with /i, /2 functions of t and x satisfying the following equation: 

— - — -/?x + 7t = 0. (25) 

If we consider the transformation ([14p between /i, /2 and the gauge function g, then equation 
([25]) becomes: 

-||-/?^ + 7i = 0, (26) 



which can be easily integrated, i.e.: 



/3 = ^(-/3x2 + 27tx), (27) 



Jacobi's original paper is in Latin, but the mathematical formulas could be understood by any mathemati- 
cian. 



Then Li can be written as follows: 



Li = i:cV + ^(-/?x2 + 27tx) + ^. (28) 



We remark that equation (j22p can be transformed into an autonomous equation by considering 
the canonical variables t, x of the only Lie point symmetry admitted for any 7 7^ 0, i.e. 

T = tdt + xdx, =^ i = log(t), 5 = |. (29) 
Then equation (j22p becomes 

dx d^x ^ o ndx /dx\ , , 

-, — —r = -px + -f-x^-3x^—-x{ — ] . (30) 

dt dt2 '^ ' dt \dtj ^ ^ 

with Jacobi last multiplier and consequently Lagrangian 

M, = eh^ =^ Li = Vx2(^y + A(t,x)^ + /2(t,x) (31) 



with /i,/2 satisfying the following equation: 

e^^Px + e^'7 - e^*x^ = 0. (32) 



dfl 9/2 ^3ir,~ , ^3i„, ^3t~3 



dt dx 

We find that the Lagrangian Li ()24p does not admit the Lie point symmetry T (j29p as a 
Noether symmetry [T3j. 

If in equation (p2|) we assume 7 = 0, i.e.: 

X X + xi + /3x = (33) 

then we find that it admits an eight-dimensional Lie symmetry algebra generated by the 



ing eight operators 
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(34) 



which means that equation (|33p is hnearizable by means of a point transformation [TT]. In 
order to find the Hnearizing transformation we have to look for a two-dimensional abelian 
intransitive subalgebra [11], and, following Lie's classification of two-dimensional algebras in 
the real plane [IT], we have to transform it into the canonical form (Type II) 

dx, idx (35) 

with i and x the new independent and dependent variables, respectively. We found that one 
such subalgebra is that generated by Fy and Fg. Then, it is easy to derive that 

i=t, x = -x^ (36) 

and equation ([33|) becomes 

d^i 

i|^0. (37) 



Now we can use the eight Lie point symmetries (|34p to generate fourteen different Jacobi last 
multipliers of equation (I33p by means of ()16p , and therefore fourteen different Lagrangians of 
equation ()33p by means of (J13p as it was shown in [22j. Here we report just three. 
The Jacobi last multiplier J35 which is derived from Tj, and Fs, i.e. 



x^ 



"^^^ {xx + pt){2txx + I3f^ - x^)' ^^^^ 



yields the Lagrangian 



^35 = —r--^—^i-(3t^-2txx + x^)log{(3t^ + 2txx-x^)-- 
2t(/3t^ -I- x^) 2i 



/3t2 + x^ 



log(/3t + xx){l3t + xx) + fi± + /2, 



with -— —- = ——-T, ^. (39) 

dt dx t{pt^ + x^) ^ ' 

The Jacobi last multiplier J58 which is derived from F5 and Fg, i.e. 

•^58 = -737^, (40) 

pt + XX 

yields the Lagrangian 



-^58 = -(/3t + xi;)log(/3t + xi;) -F/3i-Fxi;-F/ix-F/: 



2, 



wi«. f -f ^0. («) 



The Jacobi last multiplier Jyg which is derived from Tj and Tg, i.e. 

J78 = x^ (42) 

yields the Lagrangian 

L78 = ^rcV + /li; + /2, with ^-^ = px. (43) 

We note that Jjs corresponds to the multiplier ()23p found by Jacobi. 

We remind the reader that any ratio of two multipliers gives a first integral of equation ()33p . 

e.g. 

h = '^ = -x^ + 2txx + pf, (44) 

h = '■^ = -xx-l3t, (45) 

I3 = ^ = {xx + pt){x'^ -2txx-f3t^). (46) 

2.1.2 Example II 

The second example by Euler that Jacobi presented is the following equation: 

2x^x + x^i;^ - ax^ + /3t2 - 7 = (47) 

with Q, /? and 7 arbitrary constants. Jacobi obtained the multiplier 

Mi = x, (48) 

by means of (|19p and showed how to integrate equation (j47p . 

We can use Jacobi's multiplier (jl8]) to derive a Lagrangian of equation (|17|) by means of 
(UnD, i.e. 

i2 = -xi;2 + /i(t,x)x + /2(t,x) (49) 

with /i, /2 functions of t and x satisfying the following equation: 

2 f^fi 9/2 \ 2 ^.2 



If we consider the transformation ()14p between /i, /2 and the gauge function g, then: 

/3 = ^(a^' + /3t'-7), (51) 



and Lo can be written as follows: 



^2 = l""^^ + 2^^""^' + ^*' ~^^^% ^^^^ 



We note that equation (l47j) does not possess any Lie point symmetry. 
If we assume /? = 7 = 0, then equation ([^7]) . i.e. 

2x^x + x^x^ - Qx^ = (53) 

admits a two-dimensional Lie symmetry algebra generated by the following two operators 

^x = tdt + xd^, n2 = dt. (54) 

Then we can use these two symmetries to generate another Jacobi last multiplier of equation 
()53p by means of ()16p . and therefore a different Lagrangian by means of (jlSp . The Jacobi last 
multiplier J12 which is derived from f^i and ^2, i-c 



J12 = -To' (^^) 

—a + x"^ 



yields the Lagrangian 



X I X \ ( X \ 

L12 = -2^ arctanh ^ - log 1 + /li + /2, 

Va V V « / V « / 

. (9/i 5/2 1 

Finally the ratio of the two multipliers (j48p and (|55p is a first integral of equation (|53p . i.e. 

I=^^=\x{-a^±\ (57) 

If we assume a = /3 = 7 = 0, then equation (|T7|) . i.e. 

Ix^x + x^i^ = (58) 

admits an eight-dimensional Lie symmetry algebra, and fourteen different Lagrangians can be 
derived. 

3 Application of the method by Jacobi to equation ([l]) 

It is obvious that equation ([1]), for any h{x) and c(x), is a particular case of the class of 
equations (jlSp studied by Jacobi. Therefore a Jacobi Last Multiplier ()19p is already known, 
i.e. 

^^^g2P,(x)^ with Ph{x)= fb{x)dx, (59) 

10 



and the corresponding Lagrangian (j20p is 



Li = ie2^''(^)x2 + /i(t, x)x + /2(t, x) (60) 



with /i, /2 functions of t and x satisfying the following equation: 



dfl <9/2 2Ph(x) / N /ci\ 

— - — = e ^ne(,)^. (61) 



The Lagrangian derived by Museliak et al. p^, after lengthy calculations, i.e. 

L = ^e2^'(^)i;2 _ f e2P,{x)^(^)^ ^^^ (g2) 

is a subcase of the Lagrangian ([60]l . with /i = and /2 = /2(x) = f e'^^''''^'c{x)xdx, which is 
an obvious particular solution of (j6ip . Actually we can derive another Lagrangian of ([T|). In 
fact we note that ([T]) admits one trivial Lie point symmetry for any b{x) and c(x), i.e. T = dt, 
which is also a Noether's symmetry |13] for the Lagrangian (j60p . Therefore a first integral can 
be easily obtaineqj from Noether's theorem [13], i.e.: 

h = -e^^'x"^ + / e'^^''c{x)xdx. (63) 

Let us use the property (d) of the Jacobi last multiplier. If one knows a Jacobi last multiplier 
Ml in (|59p and a first integral /i in (|63p of equation ([T]) , then their product is another Jacobi 
last multiplier, i.e. 

M2 = Mih = \- e^^" (e^^^x^ + f e^^^cix)x dx\ (64) 

Consequently we are able to obtain a second Lagrangian of equation ([T]) for any b{x) and c{x), 
i.e. 

L2 = -^e'^^'x'^ (e^^^x^ + U I e'^^^c{x)xdx\ + /i(t,x)x + /2(t,x) (65) 

with /i, /2 functions of t and x satisfying the following equation: 

'c{x) I e ^c{x)xdx. (66) 



^ _ ^ - .SnW^.^N / 2P, 



dt dx 

This Lagrangian admits F = 9^ as a Noether's symmetry and the corresponding first integral 

is just the square of /i in ([631) . 

We can keep using property (d) to derive more and more Jacobi last multipliers and therefore 



^This energy-type first integral was obtained in [12] . 
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Lagrangians of equation ([T]). In fact other Jacobi last multipliers can be obtained by simply 
taking any function of the first integral /i in (j63p and multiplying it for either Mi in (j59p or 
M2 in (j64p . and so on ad libitum. For example, we may take if, and then obtain another 
Jacobi last multiplier by taking the product of Mi in (j59p and if, i.e.: 



M3 = Milf = i e^^^ (e^^^x^ + / e^^'c{x)x dxj (67) 

which yields the following third Lagrangian of equation ([T]) 

L3 = — e^^'i;^ ( e'^^'-x'^ + lOe^^'i^ / e^^''c(x)x dx + 60 I e^^''c(x)x dx j j 

+/i(t,x)x + /2(t,x) (68) 

with /i, /2 functions of t and x satisfying the following equation: 



^ _ H = e2^^(^)c(x)x n e^''^c{x)xdx\ 



(69) 



4 Final remarks 



The purpose of the present paper is to exemplify once again the method of Jacobi for finding 
Lagrangians, and to stress its strong connection with Lie symmetries and Noether symmetries. 
We advocate Jacobi Last Multiplier as an essential tool for studying nonlinear dynamical 
systems. As (apparently) stated by Henry S. Truman: 
There is nothing new in the world, except the history one does not know. 
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